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ABSTRACT 

We present our analysis of photometric data in the Johnson B and V filter of the 
southern Blazhko star SS For. In parallel, we analyzed the V observations obtained 
with the ASAS-3 photometry of the star gathered between 2000 and 2008. In the fre- 
quency spectra resulting from a Fourier analysis of our data, the triplet structure is 
detectable up to high order, both in the B and V data. Moreover, we find evidence for 
quintuplet components. We confirm from our data that the modulation components 
decrease less steeply than the harmonics of the main frequency. We derived the vari- 
ations of the Fourier parameters quantifying the light curve shape over the Blazhko 
cycle. There is good agreement between the spectroscopic abundance and the metal- 
licity determined from the Fourier parameters of the average light curve. SS For is 
peculiar as a Blazhko star because of its strong variations around minimum light. 

Key words: stars: oscillations - stars: variable: RR Lyrae - stars: individual: SS For 
- techniques: photometric. 



1 INTRODUCTION 

A large fraction of the RR Lyrae stars shows a periodic 
amplitude and/or phase modulation with a period of typ- 
ically ten to hundreds of times the pulsation period. This 
phenomenon is referred to as the Blazhko effect, after the 
Russian astronomer who first reported it (Blazhko 1907). 
The origin of the Blazhko effect is still a matter of con- 
troversy. The most widely discussed models attribute the 
effect to either the consequences of the interaction of a 
magnetic field with the main radial pulsation (Shibahashi 
& Takata 1995), or a resonance between the main mode 
and a non-radial mode of low degree (Van Hoolst, Dziem- 
bowski & Kawaler 1998; Nowakowski & Dziembowski 2001; 
Dziembowski & Mizerski 2004). However, also models that 
do not require nonradial modes have been proposed, e.g. by 
Stothers (2006). This scenario attributes the Blazhko varia- 
tion to a variable turbulent convection due to transient mag- 
netic fields in the star. All models presently proposed for the 
Blazhko effect have shortcomings in explaining the variety 
of features shown by Blazhko stars. Therefore, theoretical 
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efforts to revise or expand the models would be worthwhile, 
and even the exploration of alternative explanations. 

Since its discovery almost a century ago, most studies 
of the Blazhko effect were carried out from the northern 
hemisphere (e.g., Szeidl 1988, and references therein; Smith 
1995, and references therein). For a long time there were 
fewer well-established field Blazhko stars at southern decli- 
nations. The available data for southern field Blazhko stars 
(e.g., Hoffmeister 1956; Kinman 1961; Clube et al. 1969; 
Lub 1977) are often insufficient to determine the Blazhko 
periods with the required accuracy. Accurate and complete 
photometric data sets of southern field Blazhko stars are 
still lacking. Extensive multi-target surveys such as MA- 
CHO (Alcock et al. 2000, 2003) and OGLE (Moskalik & 
Poretti 2003), and AS AS (All Sky Automated Survey, Po- 
jmanski 2000) have significantly contributed to our knowl- 
edge of the Blazhko effect, and expanded the list of known 
Blazhko stars. Wils & Sodor (2005) published a list of new 
and confirmed Blazhko targets, and Szczygiel & Fabrycky 
(2007) obtained interesting new results, among which is a 
star with multiple Blazhko periods. 

Nevertheless, long-term campaigns dedicated to partic- 
ular Blazhko stars and yielding complete light curves at dif- 
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ferent phases in the Blazhko cycle remain of great scien- 
tific value. Only these can reveal changes in the light curve 
shape, as well as long-term changes in the characteristics of 
the Blazhko effect such as amplitudes, phases, and periods. 
In this way they can give crucial information for deciding 
among the different hypotheses for the Blazhko effect. 

Photometry with a good spread over both the pulsation 
and Blazhko cycles allows us to perform the first rigorous 
frequency analysis of a southern field Blazhko star. The first 
accurate photometric data set covering both the pulsation 
and the Blazhko cycle was published by Jurcsik et al. (2005a) 
for the northern short-period Blazhko star RR Gem. 

SS For (SAO 167572; a(J2000): 02 h 07 m 52 s , <5(J2000): 
-26° 51' 54") is a relatively bright southern Blazhko star. 
According to the General Catalogue of Variable Stars 
(Kukarkin 2003) its V brightness changes within the 
range 9.45-10.60. The General Catalogue of Variable Stars 
(GCVS) lists a period of P = 0.495432 d (Kholopov et al. 
1998) or about 11 h 53 min, corresponding to a frequency 
of f = 2.0184405 cd _1 . 

We selected this target on the basis of references to 
its variable light curve in the literature (Lub 1977). Our 
analysis of both the HIPPARCOS and the ASAS-3 photom- 
etry demonstrated the existence of a Blazhko effect with 
strong variations around minimum. The deduced Blazhko 
period of SS For (about 35 d) turned out to be short 
enough to cover several modulation cycles during one ob- 
serving season. All our data of SS For are available online 
at http://www.univie.ac.at/tops/blazhko/data/. 

In Section 2 we present the new observations as well as 
the ASAS-3 data used for comparison. The different steps in 
the data analysis and its results are described in Section 3. 
Section 4 is dedicated to a discussion of the results obtained. 
Finally, some concluding remarks are given in Section 5. 

2 OBSERVATIONS 

2.1 New photometric measurements 

Our photometric observations were carried out with 3 differ- 
ent telescopes in the southern hemisphere. All were equipped 
with photomultipliers to obtain comparative photometry in 
the Johnson B and V passbands. 

The photometric observations of SS For started from 
the South African Astronomical Observatory (SAAO) in 
Sutherland, South Africa, in October 2004, and were car- 
ried out from the 0.5-m and 0.75-m telescopes until Novem- 
ber 2005. The 0.5-m telescope was equipped with the mod- 
ular photometer. At the 0.75-m telescope we used the UCT 
photometer. From July until September 2005, the 0.6-m 
telescope at Siding Spring Observatory (SSO) near Coon- 
abarabran, Australia, joined in the observing campaign. Its 
photometer is equipped with a GaAs photomultiplier tube. 
Given the nearly 12-h period of SS For, the complementary 
longitudes of the two observatories allowed for observations 
of a larger part of the light curve. Additionally, the height 
of 1-day aliases introduced by the spectral window (see also 
Fig. 1) was decreased by the use of more than one observing 
site. This enables an easier and more accurate determination 
of the frequencies of the stellar light variations. The total 
time span of the data set is 345 days, almost 10 Blazhko 
cycles. 



Spectral Window SAAO+SSO (2004-2005) 

1 1 1 1 1 1 1 1 r 

0.9- 




0.0 1.0 2.0 3.0 4.0 5.0 

Frequency (c/d) 



Figure 1. Window Function for the data set gathered at SAAO 
and SSO (2004-2005). 



The observations were carried out using the 3-star 
technique, as described by Breger (1993). The comparison 
and check stars were HD13334 (a(J2000): 02 h 09 m 45 s , 
<5(J2000): -27° 06' 13", V =9.7, B — V = 0.57) and HD13181 
(a(J2000): 02 h 08 m 13 s , 5(J2000): -27° 01' 29", V =9.8, 
B — V = 0.51). Table 1 shows a journal of the measurements. 
Depending on the weather conditions, the accuracy of the 
individual differential measurements varied between 3 and 
10 mmag. In the reduction procedure, variable atmospheric 
extinction was taken into account. The night-to-night shifts 
of the differential magnitudes of the two comparison stars 
laid within the scatter of a few millimags. 

In order to combine the data sets from the different 
telescopes and instruments we carried out transformations 
to the standard system. This is crucial, since for variables 
with changing light curve shapes zero point offsets can of- 
ten not be correctly determined. Moreover, the large varia- 
tions in the B-V colour (more than 0.3 mag - see Figure 3) 
throughout the pulsation cycle result in a considerable dis- 
tortion of the light curves (as large as 0.02-0.03 mag) if the 
colour systems are not the same. For the SAAO data we 
used transformation coefficients determined by D. Kilkenny 
(private communication), and from our own standard star 
measurements. For the SSO data we took coefficients deter- 
mined from standard star measurements by R. Shobbrook 
(private communication) , close to the ones given by Bernd- 
nikov & Turner (2001). 

The light curves of the new data in Johnson B and 
Johnson V , folded with the main pulsation period, are shown 
in Figure 2. The differential magnitudes are shown with re- 
spect to HD13181. 



2.2 ASAS-3 data for SS For 

SS For (ASAS J020752-2651.9) was observed as part of the 
All Sky Automated Survey (ASAS-3; Pojmanski 2002) from 
November 2000 onwards. In most cases ASAS photometry 
is accurate to about 0.03 mag (Pojmanski 2000). We anal- 
ysed SS For data from ASAS-3 which were gathered be- 
tween November 2000 and July 2008 (HJD 2451868.613- 
2454656.892, 484 points over 2788 days or almost 80 Blazhko 
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Table 1. Observing log for the SS For Johnson B and V data. 



SAAO data (2004+2005) 



SSO data (2005) 
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Length [hrs] 
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N 
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SAAO 


total 


2004 
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SSO 


total (=2005) 


Observing time [hrs] 


158.32 


82.63 


75.70 


Observing time [hrs] 




43.51 




Number of datapoints 


1000 


490 


510 


Number of datapoints 




218 




Number of nights 


34 


16 


18 


Number of nights 




16 




Total time span [days] 


314 


79 


28 


Total time span [days] 




64 





Total observing time [hrs] 201.83 

Total number of datapoints 1218 

Number of nights 50 

Total time span [days] 345 



cycles). 
Fig. 4. 



The folded V ASAS-3 light curves are shown in 



ysed each of the standardized data sets (SAAO, SSO, ASAS) 
separately to check the compatibility of the frequency solu- 
tions. 



3 DATA ANALYSIS 

The frequency analyses were performed with Period04 (Lenz 
& Breger 2005), a package applying single- frequency power 
spectra and simultaneous multi-frequency sine-wave fitting. 
Fourier analyses were carried out on subsets of the photo- 
metric data in both filters. Before merging them, we anal- 



3.1 Triplet fits: may be not sufficient 

The best fits to the data were determined by means of 
a successive prewhitening strategy. From this analysis the 
triplet structures typical for RR Lyrae Blazhko stars clearly 
emerged. 

For Blazhko stars very often the a priori choice is made 
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Figure 2. B (upper panel) and V (lower panel) light curve of 
SS For, folded with the main period. The full line represents the 
mean light curve. The Blazhko effect is very pronounced around 
both minimum and maximum light. 
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Figure 3. Mean B — V color curve of SS For, resulting from 
subtracting the two mean light curves (Figure 3). 
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Figure 4. SS For data from ASAS-3, folded with the main pul- 
sation period. 




8.0 12.0 

Frequency (c/d) 

Figure 5. Fourier spectrum for the data set gathered at SAAO 
and SSO (2004-2005). The highest peak corresponds to the main 
frequency at 2.018433 cd _1 . The 1-day aliasing is clearly visible. 



to describe the light curve variations by means of equidistant 
triplets at the main frequency and its harmonics according 
to the formula (see also Kovacs 1995): 

f(t) = A + ELJ^U sin(27r/c/ t + <t> k ) 
+A k+ sin(27r(fe/ + f B )t + 4> k +) 
+A k - sin(27r(fc/ - f B )t + (j> k -)] 
+B sm(2n{f B t + ct> B )). 



(1) 



The amplitudes and phases are calculated up to order n, for 
which the amplitudes of the higher-order harmonics are still 
above the significance level. In this formula, fo is the main 
pulsation frequency and f B is the Blazhko frequency. 

In order not to presume equidistantly spaced triplets 
we also fitted the data with triplet structures without the 
condition of equidistance, i.e. the side peak frequencies are 
"let free". The resulting departure from equidistance (see 
also Breger & Kolenberg 2006) was smaller than the uncer- 
tainties on the frequencies, so we can resort to Equation (1) 
for the triplet fits. 
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We calculated optimum values for the frequencies, their 
amplitudes and phases by minimizing the residuals of the 
fit given by Equation (1). Errors on these values were deter- 
mined through extensive Monte Carlo simulations, and are 
of the same order as the error values calculated following 
Breger et al. (1999) and Montgomery & O'Donogue (1999). 
Handler et al. (2000) found that correlated noise may lead 
to an underestimation of the errors. Following their con- 
clusions we multiplied the errors resulting from our Monte 
Carlo simulations by a factor of 2, in order to get a reliable 
and realistic error margin. The uncertainties obtained and 
given in Table 2 were confirmed by our analyses of different 
subsets of the data. 

3.1.1 AS AS data 

With the ASAS-3 data set until July 2008 we find 
/o=2. 01843 ± 0.00001 cd _1 for the main frequency and f + 
/s=2.04712 ± 0.00002 cd _1 for the first detected side peak 
frequency, resulting in a Blazhko frequency / s =0. 02868 ± 
0.00002 cd _1 . Wils & Sodor (2005) published the results of a 
period analysis of a set of southern RR Lyrae stars exhibit- 
ing the Blazhko effect and published new elements for pre- 
viously known as well as new and suspected Blazhko RRab 
stars of the ASAS database. They determined the Blazhko 
periods of the variables with a fit containing the first four 
trial modulation frequencies, i.e., a fit including the triplet 
components around the main frequency fo and its harmonic 
2/o. On the basis of the earlier ASAS-3 photometry for SS 
For, they list a Blazhko period of 34.8 d, which corresponds 
well with the value we find. 

We checked the ASAS data for period changes by look- 
ing at all normal maxima in the ASAS data. 

3.1.2 SAAO and SSO data 
The V data 

The highest recorded peak-to-peak amplitude in our V 
data is 1.34 mag. The harmonics of the main radial fre- 
quency are significant up to 16/o and both triplet compo- 
nents are detected above the 3.5a noise level for combination 
frequencies up to the 11th order. 

When including the ASAS-3 V data in the analysis, we 
obtain f = 2.01843 ± 0.00001 cd _1 for the radial pulsation 
frequency and f N = 2.04708 ± 0.00009 cd _1 for the right 
side peak frequency. Given the longer time span of the data 
set, we use the frequencies found from the combined ASAS- 
3+SAAO+SSO data as start values for optimization to fit 
the new data gathered at SAAO and SSO. Note that we do 
not use the ASAS data to determine the final fit given in 
Table 2. 

The B data 

RR Lyrae stars have higher amplitudes in the blue re- 
gion of the spectrum. The highest recorded peak-to-peak 
amplitude from our B data set is 1.65 mag. For the main 
frequency and its harmonics, the amplitudes in B are on 



average a factor of about 1.23 larger than those in V for fo 
and its harmonics up to order 10. 

Since they were based upon a larger and more extended 
set of data, we used the optimized frequencies of the fit to 
the V data to fit our data obtained in the B filter. Given 
the higher amplitudes in the B filter, the harmonics of fo 
are significant up to a higher order, 17/o- Also both triplet 
components have amplitudes above the 3.5<r noise level (the 
adopted threshold for combination frequencies) up to the 
10th order. 

The final fit 

Table 2 list the results of a multifrequency fit to the 
combined set of SAAO+SSO data, according to Equation 
(1) and based on the magnitude differences relative to 
HD13181, the comparison star with coordinates closer to SS 
For and a smaller colour difference to the target. The results 
based on the magnitude differences relative to HD13334 dif- 
fer from the ones published in Table 2 within the given un- 
certainties. Note that for independent frequencies (in this 
case only fo and /b) to be significant we required a signal- 
to-noise level exceeding 4. For combination frequencies we 
required the signal-to- noise level to be higher than 3.5. 

Ephemerides 

From the new SAAO and SSO data we derived the fol- 
lowing ephemerides of maximum light and maximum pulsa- 
tion amplitude: 

HJD(T max ) = 2453296.440 + 0.495433 x £ pulsation (2) 

HJD(TBimax) = 2453298.914 + 34.94(±0.05) x Seiazhko (3) 

The time of maximum light (T max ) was recorded in our 
data. The time of the highest light curve amplitude (TBimax) 
was determined on the basis of our best fit to the data. 
Our highest recorded maximum was at HJD=2453296.440 
but according to our best fit it preceded the absolute 
maximum, which happened 5 pulsation cycles later at 
HJD=2453298.914. 

3.2 Residual power and quintuplet components 

After subtraction of the best triplet fit to the data, there is 
still residual scatter of the light curve (Figure 6, grey filled 
circles), concentrated at the phases of minumum to maxi- 
mum light (pulsation phase (f> — 0.8 — 1.1) and in the bump 
phase (pulsation phase <f> = 0.5 — 0.8). This indicates that 
there are most likely other periodicities present in the data 
set. Therefore, we searched for additional frequencies in the 
residuals after subtraction of the fit including the main fre- 
quency, its harmonics and all detectable triplet components 
(see Table 2). 

Applying this procedure to the B data in which the 
modulation components have the higher amplitudes, we first 
detected a peak at the frequency 6.1125 cd _1 with a signifi- 
cant amplitude. This peak is located at the expected position 
of the higher quintuplet component at 3/o, 3/o + 2/b. We 
also detected significant power at the following frequencies: 
12.1678 cd^ 1 (6/ + 2f B ) and 10.0349 cd" 1 (5/ - 2/ s ). 
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Table 2. Amplitudes and phases (fraction of 2tt) of the pulsation and modulation frequency components of SS For (triplets) for the best 
fit including triplet components. The values displayed in italics correspond to combination frequencies not exceeding a signal-to-noise 
level of 3.5. The residuals of the triplet fit to the V data are 0.017 mag, to the B data 0.019 mag. At the bottom of the table we list the 
most prominent quintuplet components that are found in the data set. After inclusion of the 3 prominent quintuplet components the 
residuals are 0.015 mag both in B and V . The uncertainty on /o is 10~ 6 cd _1 and the uncertainty on So + Sb is 4 X 10~ 5 cd _1 . 
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Most prominent quintuplet components found in the data: 



3/o + 2/ B 6.1125 0.008 0.54 0.04 0.014 0.85 0.03 

6/o + 2/ B 12.1678 0.008 0.88 0.04 0.011 0.97 0.04 

5/o - 2/ s 10.0349 0.008 0.26 0.03 0.009 0.33 0.03 
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Figure 6. Residuals of the B (left panel) and V (right panel) data of SS For after subtraction of the mean light curve (top panels, open 
circles) and after subtraction of the best fit including triplets (top and bottom panels, filled grey circles). When including the quintuplet 
components in the fit (bottom panels with a larger magnitude scale, filled black circles), the residuals reduce even more. 
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20.0 



Figure 7. Fourier spectrum showing the quintuplet components 
found in the data for SS For after subtraction of the main fre- 
quency, its harmonics and the triplet components. The notation 
"9fe+" and stands for fc/o + 2/s and fe/n — 2/g respectively; 

"a" denotes alias peaks. 



separations within the uncertainty on the frequency. As 
it is quite remarkable to find additional frequencies exactly 
at the expected positions of quintuplet components, we have 
added the found frequencies, their amplitudes and phases 





12.12 12.16 12.20 12.24 

Frequency (c/d) 



Figure 8. Fourier spectrum for an additional frequency compo- 
nent found in the B data after subtraction of the triplet fit plus 
the found quintuplet components. Left panel: Spectrum between 
0-20 cd _1 , right panel: zoom. The peak is at the expected position 
of a scptuplct component (6/o + 3/a). 



below Table 2. Figure 7 shows the Fourier spectrum in the 
B data after subtraction of the triplet solution. Quintuplet 
components and aliases are indicated. We clearly see even 
more quintuplet components in the Fourier spectrum than 
were obtained through subsequent prewhitening. The quin- 
tuplet components emerge most clearly from the B data set, 
but we find the same frequencies with significant amplitudes 
in the V data. 

After including the detected quintuplet components in 
the fit, we find residual power at several frequencies around 
the significance level. None of them is clearly a quintuplet 
component and at this level spurious peaks occur. More 
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quintuplet components may be present in the data, but were 
not unambiguously detected by us. Given the imperfect cov- 
erage of our data (gaps) we considered it inappropriate to fit 
more quintuplet components to the data. We note, however, 
that in the B data we find evidence for a significant peak 
at f2.f964 cd _1 , corresponding to (6/0 + 3/s, see Figure8). 
More evidence is needed to confirm this frequency as one re- 
ally present in the data and not an artefact of the sampling. 
Jurcsik et al. (2008) identified frequencies at kfo ± 4/s in 
data of MW Lyr. They also found septuplet frequencies (i.e., 
frequencies at kfo ± 3/s) and kfo ± 5/s components in the 
residuals of their data set, though much less significantly. 
Multiplet components kfo ± j'/s with j ' ^ 3 (i.e., septuplet, 
nonuplet, etc. components) have not been considered in any 
of the current models for explaining the Blazhko effect. 

Table 2 gives the triplet components, as these produce 
a stable fit and hence are suited for calculations of derived 
quantities (with appropriate error bars). We also applied a 
fit including equidistant quintuplet components to the data 
set, but found this fit to be unstable. For the main frequency 
components (including harmonics) there is no significant dif- 
ference between the triplet and quintuplet fit. In the fit in- 
cluding all quintuplet components, however, some side peak 
amplitudes and phases may drift because of the high number 
of free parameters and close frequencies. 

As shown in Figure 6, the scatter in the residuals de- 
creases when the quintuplet frequencies are taken into ac- 
count. The new fit including the detected quintuplet compo- 
nents makes it possible to more accurately describe the vari- 
ations at the bump phase and around minimum and maxi- 
mum light. 

3.3 Properties of the modulation components 

The properties of the modulation components kfo + /s and 
kfo — fs in the frequency spectra of a Blazhko star consti- 
tute an important test for the models proposed to explain 
the modulation. However, it is important to obtain reliable 
uncertainties for the amplitudes and phases of the modula- 
tion components. As was shown by Jurcsik et al. (2005a), 
the amplitudes of the triplet (or generally: multiplet) com- 
ponents depend strongly on the coverage. 

3.3.1 Amplitude ratios, phase differences and asymmetry 
parameters 

From their extensive study comprising 731 Blazhko variables 
in the LMC, Alcock et al. (2003) found that the relative 
amplitudes of the first order modulation components are 
usually in the range 0.1 < A±/A\ < 0.3, but there are 
some extreme values. For our V data, we find A + /A\ — 
0.118 ± 0.007 and A- j A x = 0.057 ± 0.007. For the B data, 
A + /A 1 = 0.111 ± 0.005 and = 0.053 ± 0.005. 

Table 3 lists the amplitude ratios, phase differences and 
their errors for the first six modulation component pairs, 
defined as 

Rk = A k f 0+ f B /A kfo -f B (4) 
and 

A<j) k = 4>kf +f B - <t>kf -f B ■ (5) 



We also list the so-called asymmetry parameter 



introduced by Alcock et al. (2003) to quantify the degree 
of asymmetry in the peaks. The distribution of the Q pa- 
rameter for the Blazhko stars from the MACHO data base 
(Alcock et al. 2003) peaks at +0.3. 

According to Kovacs (1995, and references therein), the 
parameter combinations R k and A<j)k are expected to be 
constant under the assumption of a simple oblique rotator- 
pulsator model. However, such a simple model seems inad- 
equate to explain the variety of observations related to the 
Blazhko effect, and hence we give Table 3 rather as light 
curve diagnostics (see also Smith et al. 1999 and Jurcsik et 
al. 2005a). 

Due to the rather large error bars on the parameter 
combinations, the phase differences A(/> k , especially for the 
V data, do not deviate significantly from a constant value. 
On the other hand, the amplitude ratios R k show large 
changes, both in V and in B, and with the obtained error 
bars they can definitely not be considered as constant. 

As can be seen in Table 3 the Q value can change its 
sign in different orders. The positive Q values at multiplet 
order k = 1, 4 and 5 point to an asymmetry with higher 
amplitudes at the higher frequency lobes in the triplets, as 
is mostly the case in Blazhko stars. At multiplet orders k — 
2 and 3 the left side peaks are higher, yielding negative Q 
values. 



3.3.2 Modulation in B versus V 

The Blazhko modulation in the B filter is larger than in the 
V filter. Jurcsik et al. (2005b) found, from a combination 
of different data sets, that the value of J 4 mod (i3)/A mod (l / ) 
lies within the range of 1.23 - 1.39 with a mean value 
of 1.30. Their value was determined from the sum of the 
Fourier amplitudes of the first four modulation components, 

A fo+fs > A fo-f B > A 2f +f B > A 2f -f B ■ For our data we obtain 
a value A ulod (B)/A rno<1 (V) = 1.22 ± 0.10, taking into ac- 
count the first four modulation components. At 2/o + /s 
the side peak amplitude determined from our data set is 
rather low. When the higher order side peaks are taken into 
account the value of A nlo a(B)/A lno< i(V) does not increase 
significantly (see Table 4). This implies that the modulation 
is indeed stronger in B than in V, and by the same factor 
as holds for the main frequency and its harmonics (Section 
3.2.1). 

3.3.3 Amplitude decrease for the side peaks 

Jurcsik et al. (2005a) reported that the decrease in am- 
plitude of the successive harmonics of fo is much steeper 
than that over the side peaks. Whereas the latter decrease 
can be described as exponential, the side peak amplitudes 
have a more linear decrease. As a consequence, the relative 
contribution of the modulation at higher orders is larger. 
This is illustrated in Figures 9 and 10, showing the am- 
plitude ratios A(f k )/A(f Q ), A{f k + f B )/A(f + f B ), and 
A (fk — /b)/-A(/o — fs) for the V and the B data respec- 
tively. The amplitudes of the left and right side peaks de- 
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Table 3. Side lobe amplitude ratios R k , phase differences A(p k (fraction of 2tt) and asymmetry parameters Q as defined in the text, and 
their respective errors for the V and B data, k denotes the multiplet order. 



k 


R k (V) 




A^ fe (V) 




Q(V) 




R k {B) 








Q(B) 


°Q(B) 


1 


2.05 


0.28 


0.15 


0.02 


0.34 


0.07 


2.09 


0.24 


0.14 


0.02 


0.35 


0.06 


2 


0.43 


0.24 


0.18 


0.07 


-0.40 


0.31 


0.26 


0.16 


0.25 


0.09 


-0.58 


0.50 


3 


0.53 


0.19 


0.19 


0.04 


-0.30 


0.16 


0.54 


0.17 


0.17 


0.04 


-0.30 


0.14 


4 


1.69 


0.51 


0.12 


0.04 


0.26 


0.11 


1.68 


0.39 


0.12 


0.03 


0.25 


0.08 


5 


1.18 


0.36 


0.11 


0.04 


0.08 


0.04 


1.08 


0.27 


0.09 


0.03 


0.04 


0.01 


6 


1.13 


0.46 


0.13 


0.05 


0.06 


0.03 


1.00 


0.32 


0.12 


0.04 


0.00 


0.00 



Table 4. The ratio of the strength of the modulation components 
in the Johnson B and V filter for increasing order k. k = 1 takes 
into account the side peak amplitudes around /n,...,fc/n. 



Order k 


Anod(-B) 


A mod (V) 


A mod (B)/A mod (y) 


1 


0.091 ± 0.004 


0.075 ± 0.004 


1.21 ± 0.09 


2 


0.115 ± 0.006 


0.094 ± 0.006 


1.22 ± 0.10 


3 


0.146 ± 0.007 


0.121 ± 0.007 


1.20 ± 0.09 


4 


0.185 ± 0.008 


0.152 ± 0.008 


1.21 ± 0.09 


5 


0.219 ± 0.009 


0.180 ± 0.009 


1.22 ± 0.08 


6 


0.246 ± 0.010 


0.201 ± 0.010 


1.22 ± 0.08 


7 


0.268 ± 0.011 


0.219 ± 0.011 


1.22 ± 0.08 


8 


0.282 ± 0.012 


0.231 ± 0.012 


1.22 ± 0.08 


9 


0.294 ± 0.013 


0.241 ± 0.013 


1.22 ± 0.08 


10 


0.304 ± 0.013 


0.248 ± 0.013 


1.23 ± 0.09 
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Figure 10. Amplitude ratios A(f k )/A(f ), A(f k + f B )/A(f + 
/b), and A(f k — /b)/A(/o — fs) of the detected frequencies for 
the B data. See also caption of Figure 9. 



curve (open circles). It shows that the Blazhko modulation 
in SS For is not only concentrated in the minimum to max- 
imum phase interval, as is the case for other Blazhko stars 
(see, e.g., Jurcsik 2005a). The strong variations in the bump 
phase before minimum light (<p — 0.5 — 0.8, see also Guggen- 
berger & Kolenberg 2006) are reflected in the scatter. As is 
known for Blazhko variables, the scatter is generally largest 
around the phase of maximum light. For SS For, there is a 
comparable scatter around minimum light. A large modula- 
tion around minimum light was also detected in the Blazhko 
star MW Lyr (Figure 3 in Jurcsik et al. 2008). 



Figure 9. Amplitude ratios A(f k )/A(f ), A(f k +f B )/A(f +f B ), 
and A(f k — /s)/A(/o — /s) of the detected frequencies for the V 
data. The amplitudes of the main frequency and its harmonics de- 
crease exponentially with increasing order. The amplitudes of the 
left and right side peaks, on the other hand, decrease less steeply 
and show more irregular behavior at lower multiplet orders. 

crease less steeply. They also show more irregular behavior 
at lower multiplet orders (k = 2,3,4). 

3.3-4 Light curve modulation 

Figure 6 shows the residual light curve of SS For folded with 
the main pulsation period after subtraction of the mean light 



3.4 Variations over the Blazhko cycle 

In order to assess the light curve variations over the Blazhko 
cycle, our data were divided into 10 subsets. Each subset 
corresponds to a 0.1 phase interval over the Blazhko cycle. 
For the ephemeris of Blazhko phase $ = 0, defined as the 
phase of the brightest maximum, we adopted the time of the 
highest detected maximum obtained from the best fit to our 
observations, i.e., T = 2453298.914. 

Table 5 gives the time base for each of these intervals, as 
well as the number of data points. For the B data the gaps 
in the light curves at different Blazhko phases do not allow 
us to perform a fit for each of the phase intervals (Figure 14, 
available online) . If coverage of minimum or maximum light 
is lacking, it is impossible to reliably fit the data and de- 
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Table 5. A log of the subsets of the observations used to construct the light curves at different phase intervals in the Blazhko cycle. For 
the V data set we used SAAO, SSO and ASAS data, for the B filter only SAAO and SSO. 



A* 


HJD (-2450000) (V) 


N(V) 


HJD (-2450000) (B) 


N(B) 


0.0-0.1 


1868.61332 


- 4629.92251 


201 


3296.440-3612.133 


132 


0.1-0.2 


1870.60972 


- 4527.50794 


248 


3300.380-3581.324 


184 


0.2-0.3 


1873.60768 


- 4633.91429 


85 


3303.413-3584.678 


174 


0.3-0.4 


1878.59277 


- 4534.49935 


84 


3587.590-3622.167 


38 


0.4-0.5 


1880.57795 


- 4501.57064 


133 


3590.204-3591.321 


24 


0.5-0.6 


1884.56548 


- 4645.90353 


163 


3592.519-3594.678 


96 


0.6-0.7 


1888.56074 


- 4649.89091 


241 


3319.402-3633.294 


137 


0.7-0.8 


2067.93681 


- 4652.92994 


209 


3288.278-3602.671 


204 


0.8-0.9 


1930.5492 


- 4656.89163 


110 


3289.481-3571.323 


158 


0.9-1.0 


1900.53501 


- 4520.53379 


206 


3296.257-3574.686 


66 



termine Fourier parameters. For the V data, however, we 
added the ASAS-3 photometry, and fitted each of the phase 
intervals with a higher-order harmonic fit according to: 



f(t) = A + ^ A k sin(27rfc/ fc t + fa), 



(7) 



and calculated Fourier parameters, namely the amplitude 
ratios 



Rki = A k /A 1 

and the epoch-independent phase differences 

4>kl = 4>k — kfa. 



(8) 



(9) 



These parameters offer a way to quantify the shape of the 
lightcurves. Fig. 12 shows the variations of the amplitudes 
Ak and the phases <j)k for k = 1, 4, and Figure 13 the de- 
rived Fourier parameters. The Fourier parameters are also 
listed in Table 6. It is clear that the first order amplitude Ai 
decreases towards minimum Blazhko amplitude more dras- 
tically than do the higher order amplitudes A2, A3, A4. As 
a consequence, the Rki parameters reach their maximum 
amplitude around minimum light. Small variations in the 
fa, and subsequently the fal parameters, can be observed 
(Figure 12). The variation of fa is the largest, as was also 
noted by Jurcsik et al. (2008). 

Our observed Fourier parameters for SS For clearly fall 
within the ranges of Fourier parameters derived from the 
V light curves of 257 RRab stars published by Kovacs & 
Kanbur (1998). The ranges we obtain also intersect with the 
theoretical parameters calculated by Feuchtinger (1999) and 
Dorfi & Feuchtinger (1999, their Figure 8 and 11). Note that 
the latter authors express their Fourier parameters in the 
cosinus frame. For the amplitude ratios we have values in the 
range R 21 = 0.41-0.55, R31 = 0.27-0.34, R 41 = 0.15-0.25 
for the V data. The phase differences in the cosinus frame 
are in the ranges fax = 3.67 - 4.06, fai = 1.44 - 1.99, 
04i = -0.82 - 0.05. 
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Figure 12. Variations of the amplitudes A k and the phases <j>k 
over the Blazhko cycle for k = 1, ...,4. 



4 DISCUSSION 
4.1 Blazhko period 

The Blazhko character of SS For was recognized by Lub 
(1977), who noticed strong variations around minimum light 
of the star. At that time, its Blazhko period was not known. 



Wils & Sodor (2005) published a study on Blazhko vari- 
ability from the ASAS database. SS For appears in their 
table with elements of new and suspected Blazhko RRab 
stars. They determined a Blazhko period of 34.8 d. From 
our new data set we find a value for the Blazhko period 
Pb — 34.94 ± 0.05 d, as the difference between fo + /b and 
/o- 
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Figure 11. V light data (crosses) in different 0.1 phase intervals of the Blazhko cycle (34.94 ± 0.05 d from this data set). We used the 
SAAO, SSO and ASAS data for constructing the phase bins. The mean light curve over all Blazhko phase intervals is also shown as a 
solid line. Note the strong variations at the bump phase before minimum light. 



4.2 Multiplets challenge the models 

The observed asymmetry of the modulation components in 
the triplets, also observed in SS For (see Section 3.3.1), re- 
mains unexplained by both the resonance model and the 
magnetic model, at least in the degree it is observed. 

In our data we find clear evidence for quintuplet com- 
ponents in the vicinity of the harmonics of the main fre- 
quency in our data. A quintuplet structure is generally pre- 
dicted by the magnetic model (Shibahashi & Takata 1995), 
where the magnetic field causes the main radial mode to de- 
form and have additional quadrupole (I — 2) components. 
For more than a decade after the magnetic model was first 



presented, no evidence for quintuplet components was found 
from any data set. The argument given by proponents of the 
magnetic model was that a quintuplet structure may man- 
ifest itself as only a triplet depending on the geometrical 
configuration (angles of pulsation axis, magnetic axis and 
aspect angle). Also, the quintuplet components might have 
such a low amplitude that they are hidden in the noise of 
the frequency spectrum. This is supported by recent findings 
by Hurta et al. (2008) in RV UMa, Jurcsik et al. (2008) in 
MW Lyr, and this work, but it does not imply the magnetic 
model is the one to be preferred over the resonance model. 

A magnetic field of about 1 kG is needed in the mag- 
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Table 6. Variations of the Fourier parameters 


M, Rki , 


01, and 4>kl ovcr the Blazhko cycle for k = 1, 4. 




A* Ax R21 


R31 


-R41 01 021 031 


041 



0.0-0.1 0.502 ± 0.006 0.43 ± 0.02 0.27 ± 0.01 0.19 ± 0.01 -0.05 ± 0.01 2.28 ± 0.05 4.78 ± 0.07 1.07 ± 0.09 

0.1-0.2 0.480 ± 0.006 0.42 ± 0.01 0.28 ± 0.01 0.20 ± 0.01 -0.01 ± 0.01 2.18 ± 0.04 4.64 ± 0.06 0.91 ± 0.08 

0.2-0.3 0.450 ± 0.017 0.43 ± 0.03 0.29 ± 0.03 0.20 ± 0.03 0.06 ± 0.03 2.10 ± 0.10 4.58 ± 0.13 0.75 ± 0.20 

0.3-0.4 0.411 ± 0.009 0.41 ± 0.03 0.29 ± 0.03 0.22 ± 0.03 0.07 ± 0.03 2.11 ± 0.08 4.60 ± 0.12 0.92 ± 0.16 

0.4-0.5 0.340 ± 0.008 0.46 ± 0.04 0.34 ± 0.05 0.25 ± 0.03 0.04 ± 0.04 2.23 ± 0.11 5.06 ± 0.14 1.42 ± 0.22 

0.5-0.6 0.352 ± 0.007 0.55 ± 0.02 0.34 ± 0.02 0.17 ± 0.02 0.04 ± 0.02 2.26 ± 0.06 4.86 ± 0.08 1.26 ± 0.14 

0.6-0.7 0.374 ± 0.005 0.54 ± 0.02 0.30 ± 0.02 0.15 ± 0.02 -0.04 ± 0.02 2.36 ± 0.06 4.93 ± 0.08 1.58 ± 0.14 

0.7-0.8 0.402 ± 0.009 0.50 ± 0.02 0.33 ± 0.02 0.16 ± 0.02 -0.09 ± 0.02 2.49 ± 0.06 5.13 ± 0.08 1.62 ± 0.15 

0.8-0.9 0.458 ± 0.010 0.46 ± 0.02 0.31 ± 0.02 0.18 ± 0.03 -0.07 ± 0.02 2.38 ± 0.06 5.01 ± 0.09 1.36 ± 0.15 

0.9-1.0 0.485 ± 0.005 0.45 ± 0.01 0.32 ± 0.01 0.21 ± 0.01 -0.08 ± 0.01 2.33 ± 0.03 4.95 ± 0.05 1.20 ± 0.07 
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Figure 13. Variations of the Fourier parameters R^i and 0j,i 
over the Blazhko cycle for k = 1, ...,4. 



netic model for the amplitude modulation to be observ- 
able (Shibahashi & Takata 1995). The only Blazhko star for 
which a dedicated spectropolarimetric campaign has been 
carried out so far is RR Lyr, the prototype of the class, and 
the results are rather contradictory. Babcock (1958) and Ro- 
manov et al. (1994) reported a variable magnetic field in RR 
Lyr with a strength up to 1.5 kG, whereas Preston (1967) 
and Chadid et al. (2004) contradict these measurements. 
The detection of magnetic fields in RR Lyrae stars has been 



hampered by the fact that these stars are quite faint; RR 
Lyr is by far the brightest with V = 7.2-8.2. With the new 
generation of powerful spectrographs attached to big tele- 
scopes, however, many more RR Lyrae stars have now come 
within reach for spectropolarimetry. 

Interestingly, in this data set we only find clear evi- 
dence for quintuplet components around the harmonics of 
the main frequency and not around the main frequency it- 
self. This may be due to the fact that the observed am- 
plitude spectrum is also influenced by the coverage of our 
measurements (see Jurcsik et al. 2005a). Satellite data with 
quasi-uninterruped coverage over several Blazhko cycles will 
give an answer to many questions concerning the frequency 
spectrum of Blazhko stars, the modulation components and 
the occurrence of quintuplet components. 



4.3 Physical parameters from the light curve 

Jurcsik & Kovacs (1996), herafter JK96, determined a 
P ~ c/»3i — [Fe/H] relation for fundamental mode RR Lyrae 
stars. Their calibration was based on a sample of 81 RRab 
field Blazhko stars, and its application to RR Lyrae stars 
with independent spectroscopic metallicities has proven that 
it has an overall prediction accuracy of 0.12 dex (Kovacs 
2005). The question, however, is whether this empirical for- 
mula is also applicable to amplitude- and phase-modulated 
fundamental mode RR Lyrae stars. 

According to Layden (1994), SS For has a spectroscopic 
iron abundance [Fe/H] spcc = -1.35, which is -1.09 on the scale 
of the JK96 relation. Based on ASAS data, Kovacs (2005) 
calculated the iron abundance [Fe/H]F OU r according to the 
JK96 formula, and found a A[Fe/H] = [Fe/H] sp ec-[Fe/H] Four 
of 0.44. Kovacs (2005) attributed the large deviation both 
to an inaccurate spectroscopic iron abundance determina- 
tion, as well as the light curve shape and therefore Fourier 
parameter variations due to the Blazhko effect. We support 
the latter explanation, since from our determination of the 
Fourier parameters at different phases in the Blazhko cycle, 
we obtain values for 03i between 4.58 and 5.13 rad (in the 
sinus frame, as required for the JK96 formula). This leads 
to metallicities [Fe/H]p our between -1.55 and -0.81, a rather 
large range but with an average (-1.18) rather close to the 
spectroscopic value. 

Based on the result of Alcock et al. (2003), it is pos- 
sible that the average light curves of Blazhko stars can be 
employed equally well in the empirical formula calibrated 
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on strictly periodic stars. However, long-term monitoring - 
from months to years - is needed to achieve sufficient cov- 
erage for Blazhko variables. The test can be done for SS 
For, for which we obtained average light curve parameters 
with a high accuracy. For the mean light curve, this yields 
the value ^31 = 4.935 ± 0.057 (this corresponds to 0.785 ± 
0.009 in fractions of 2n, as can be derived from Table 2). Ap- 
plication of the JK96 calibration then yields [Fe/H]=-1.07 
± 0.08, agreeing very well with the value -1.09 by Layden 
(1994) on the JK96 scale. 

4.4 Variations at minimum light 

While in most Blazhko stars the variations of the pulsa- 
tion maximum are most pronounced, in SS For also the 
minimum shows strong variations during the Blazhko cycle. 
These variations have already been reported by Lub (1977) 
and have been investigated in more detail by Guggenberger 
& Kolcnbcrg (2006) on the basis of the data presented in 
this paper. They found a periodic behaviour of the so-called 
bump in the light curve that appears just before minimum 
light. The bump occurs at an earlier phase around Blazhko 
minimum and later around Blazhko maximum. Also the 
strength of the bump shows a dependence on the Blazhko 
phase: the bump is most distinct during minimum and al- 
most vanishes during Blazhko maximum. 

The bump has been explained by a collision between 
the free-falling high atmospheric layers and the deep atmo- 
sphere, which has a smaller infall velocity (Gillet & Crowe 
1988). From spectroscopic data of RR Lyr, Preston, Smak, & 
Paczynski (1965) concluded that a displacement of the shock 
forming region occurs over the Blazhko cycle. This may in- 
duce an earlier or later occurence of the bump in the light 
curves. Also the U-B excess observed on the rising branch, 
as well as the strength of the H emission (Hardie 1955), both 
linked to the shocks in RR Lyrae stars, are variable over the 
Blazhko cycle (Preston, Smak, & Paczynskil965). 



5 CONCLUSIONS 

High-precision photometric observations in the Johnson B 
and V filter were obtained from the South African Astro- 
nomical Observatory (SAAO) and Siding Spring Observa- 
tory (SSO), Australia, between October 2004 and September 
2005. The time span of the data covers almost 10 complete 
Blazhko cycles. 

A detailed analysis of the combined standardized pho- 
tometric measurements led to the following conclusions: 

• From a Fourier analysis of the combined SAAO+SSO 
data set we obtain f = 2.018433 ± 0.000001 cd" 1 for the 
radial pulsation frequency and /at = 2.04705 ±0.00004 cd _1 
for the right-side peak frequency. 

• From our data we find a value for the Blazhko period 
P B = 34.94 ±0.05 d. 

• The amplitudes in the B data are a factor of about 
1.2 larger than in the V data. The harmonics of the main 
frequency as well as the triplet structure around the main 
frequency and its harmonics are detectable up to high or- 
der, especially in the B data (15th order). Also the Blazhko 



modulation is a factor of 1.2 larger in the B than in the V 
data. 

• After subtracting the fit including triplet components, 
we find clear evidence of quintuplet frequencies in the 
data. We found several clearly significant peaks located at 
fc/ ±2/s frequency values. These frequencies emerged most 
clearly from the B data set in which the modulation compo- 
nents have higher amplitudes, but were also retrieved from 
the V data. Our data thus clearly show the existence of quin- 
tuplet components in SS For, as were also found by Hurta 
ct al. (2008) in RV UMa and by Jurcsik ct al. (2008) in 
MW Lyr. Quintuplet components are predicted by the mag- 
netic model for explaining the Blazhko effect (Shibahashi 
& Takata 1995), but so far there has been no convincing, 
unambiguous positive detection of a magnetic field in RR 
Lyrae stars. 

• We also find evidence for a so-called septuplet compo- 
nent in our data (at 6/o + 3/s). Multiplet components of 
order j ^ 3 in kfo ± jfs have also been found by Jurcsik et 
al. (2008) in their extensive data set of MW Lyr. They have 
yet to be explained by the models for the Blazhko effect. Un- 
interrupted satellite data of Blazhko stars, covering several 
Blazhko cycles, such as delivered by the COROT satellite 
(Baglin 2007), will undoubtedly shed new light upon the ex- 
istence of higher-order multiplet structures in Blazhko stars. 
We may expect to find many more in high-quality data sets. 

• A subdivision of the data into 10 Blazhko phase in- 
tervals shows the variations of the light curve over the 
Blazhko cycle. As all the Blazhko phase intervals in V have 
sufficient coverage, we calculated the Fourier parameters 
and plotted their variations. Observed values at different 
phases of the Blazhko cycle may be useful for confronta- 
tion with non-linear convective models such as those devel- 
oped by Feuchtinger (1999). It would be worthwhile to find 
out whether the observed variations in the Fourier parame- 
ters can be reproduced theoretically for a star with constant 
mass and metallicity. 

• Application of the empirical P — <f>:n~ [Fe/H] relation 
developed by Jurcsik & Kovacs yields a good result for the 
average light curve derived from all our data. This strength- 
ens the assumption that the average light curve of Blazhko 
stars can be employed in the empirical formula calibrated 
by Jurcsik & Kovacs on stricly periodic stars. However, the 
average light curve of a Blazhko star can only be determined 
accurately if the Blazhko cycle is sufficiently covered by the 
available data. 

• From a study of the light curve shape at different 
Blazhko phases we clearly see that the strong variations 
around minimum light in SS For are related to the position 
and strength of the bump. A detailed study of the bump 
behaviour over the Blazhko cycle, in photometric as well 
as spectroscopic data, may shed new light upon the under- 
standing of the Blazhko effect. 

Theoretical efforts to revise or expand the existing models 
for the Blazhko effect would be worthwhile, as well as the 
exploration of alternative explanations. 
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Figure 14. B light data (crosses) in different 0.1 phase intervals of the Blazhko cycle (34.94 ± 0.05 d from this data set). The mean 
light curve over all Blazhko phase intervals is also shown as a solid line. Note the strong variations at the bump phase before minimum 
light. 
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